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The vacuum stability bounds on flavour-violating trilinear soft terms are revisited from the view- 
point that one should not ban a standard-model-like false vacuum as long as it is long-lived on 
a cosmological timescale. The vacuum transition rate is evaluated numerically by searching for 
the bounce configuration. Like stability, a metastability bound does not decouple even if sfermion 
masses grow. Apart from being more generous than stability, the new bounds are largely inde- 
pendent of Yukawa couplings except for the stop trilinears. With vacuum longevity imposed on 
otherwise arbitrary LR insertions, it is found that a super flavour factory has the potential to probe 
sparticle masses up to a few TeV through B and r physics whereas the MEG experiment might 
cover a far wider range. In the stop sector, metastability is more restrictive than any existing ex- 
perimental constraint such as from electroweak precision data. Also discussed are dependency on 
other parameters and reliability under radiative corrections. 



I. INTRODUCTION 

Flavour physics is an important means to test any the- 
oretical description of elementary particles. Among nu- 
merous proposals, the one still most compelling would be 
the Minimal Supersymmetric Standard Model (MSSM) 
that offers solutions to problems not addressed by the 
Standard Model (SM). In the MSSM, the Lagrangian 
has soft supersymmetry breaking terms. Among them, 
the scalar mass and the trilinear terms are two rep- 
resentative sources of flavour/CP violation apart from 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Un- 
less tuned in a particular way, they could add to an ex- 
isting flavour-changing-neutral-current (FCNC) process 
such as B — > X s "f or create a new one, e.g. /i — > ej, 
through loop corrections. In this way, flavour physics 
provides an invaluable input to the construction of an 
underlying theory that implements spontaneous super- 
or flavour symmetry breaking. 

For the same reason, flavour physics has been a tradi- 
tional tool to unveil supersymmetry. Now that the Large 
Hadron Collider (LHC) has started its successful opera- 
tion, the significance of this role has diminished. Nev- 
ertheless, there are at least two areas in which flavour 
studies could make contributions. The first is to scruti- 
nise the flavour structure of newly discovered particles. 
This is a complement to the LHC that is less sensitive 
to extra flavour /CP violation. The second is to continue 
the search for supersymmetric particles. This should be 
vital particularly in a scenario where sparticle masses are 
around the LHC reach or beyond. In either case, the 
crucial quantity is the size of a flavour- violating term in 
the Lagrangian, often represented by the dimensionless 
quantity S to be defined in the next section. 

According to the chirality structure, a mass insertion 
(MI) can be classified as either chirality-preserving or 
chirality-flipping. In a high sparticle mass region, the 
latter type becomes generically more important than 
the former pQ. Taking for example the ratio of the 
one-loop B — >■ X s ^ amplitudes caused by the LR and 
the LL insertions, one finds a quantity proportional to 



(M 3 /m b ){(S% 3 ) LR /(Sl 3 ) LL } (see e.g. Ref. 0). It shows 
that the relative contribution from the LR insertion 
grows linearly as the gluino mass M 3 . This is one of 
the reasons why the subject of this study is the trilinear 
terms. 

The obvious problem is that one needs to get an idea of 
the flavour structure of particles that might yet be found. 
A popular strategy is to employ a model or ansatz that 
leads to a particular pattern of the soft terms. There 
are many works in this approach, and this article is not 
going to be another. The plan here is the following. The 
most optimistic scenario is presumed in which any LR 
insertion can take an arbitrary value as long as it stays 
within existing experimental and theoretical limits. Us- 
ing this input, one can estimate the maximum possible 
new signal in an FCNC process. Then it is compared 
with sensitivities of running or forthcoming experimen- 
tal searches. 

As for the input mentioned above, one would learn 
nothing new if only the present FCNC limits were used. 
One needs information from outside flavour physics. 
There are such theoretical bounds that arise from vac- 
uum stability. If a trilinear term is too large, the 
MSSM scalar potential develops a charge-and/or-colour- 
breaking (CCB) vacuum deeper than the the standard- 
model-like (SML) local minimum 014] or an unbounded- 
from-below (UFB) direction in the field space Most 
notably, the conceptual design report of a super flavour 
factory (SFF) has made an assessment of its ability to re- 
construct a LR squark insertion, imposing these bounds 
0. 

However, it is questionable whether the vacuum stabil- 
ity is an essential requirement in particular when one is 
performing an analysis of generic soft terms that is sup- 
posed to be objective. More sensibly, one could use the 
lifetime of the SML vacuum as the criterion. This view- 
point would not make sense if the universe were driven 
away from there in its history. In the case of a squark 
direction, its temperature-dependent mass [6 may defer 
the formation of a CCB minimum until the electroweak 
phase transition has happened [3 [8]. Supersymmetry 
breaking by inflation provides a more general mechanism 



2 



to lift flat directions so that relaxation to the symmetric 
phase can take place [9]. The bottom line is that there 
seems to be no definite reason to preclude a lower energy 
CCB state just for its existence. In this approach, there 
have been attempts to obtain a viable parameter space 
of the flavour-conserving stop trilinear [10] or the stau 
trilinear coupling [TT] , the sparticle [HI [13] or the Higgs 
soft masses [14] . or messenger couplings for gauge medi- 
ation [15] . This work is an application of the same idea 
to the flavour-violating trilinear couplings. 

The rest of the article begins with an abridged account 



of theoretical backgrounds in section [TTl In section III 



more practical details of the analysis method are given. 
The metastability bounds thus obtained are presented in 
section |IV| in which their behaviours are also explained. 
Section |V| is devoted to an application of the findings 
to flavour physics. Finally, the article is summarised in 
section [VT] 



II. THEORY 

For concise presentation, the following notations shall 
be employed. The species of a given matter field is de- 
noted by / which can be one of e, d, it, standing for 
charged leptons, down-, and up-type quarks, respec- 
tively. A field symbol without the superscript c refers 
to a component of the SU(2) doublets, L — (i/,e),Q = 
(u,d),H d = (H°,Hj ),H u = (H+,H°). Also, another 
redundant symbol for the down-type Higgs doublet is in- 
troduced, i.e. H e = H<i- 

The scalar potential of the MSSM is in the form 



V = V D + V F + V S 



soft 1 



(1) 



where the the F-, and the soft term contributions 
read 



V f = 



E 



dW 



9<t> a 



(2) 
(3) 



V soft = L}(Ml)ij Lj + e?(M*.)ij e 



+ QliM&ij Qj + df(M% a )ij d) + uf(M^a)ij u) 
+ 2 Re [H d L t A e ZJ e) + H d Q t A% d) - H u Q t A" S£] 
i 2 Hu \H u \ 2 + 2Re[bH u H d ], (4) 



m 2 Hd \H d \ 2 



For discussion of supersymmetric flavour violation, it 
is convenient to use the language of mass insertion ap- 
proximation [IB]. In order to specify the families of a 
given species /, one should take the super-CKM basis 
that leads to the diagonal Yukawa matrix, \{j — Xf^Sij 
with positive eigenvalues. In this article, a LR mass in- 
sertion is related with a trilinear coupling by 
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Note that this definition of average mass is different from 
that in Ref. [2], i.e. the geometric mean. Nevertheless, 
this variance does not matter since the soft scalar masses 
shall be assumed to be nearly degenerate. A trivial but 
useful identity obeyed by the LR insertions is (${j)LR — 

(dji)*RL' which follows from the Hermiticity of a sfermion 
mass matrix. 

It should be instructive to review the stability bounds 
on flavour- violating trilinear couplings as the poten- 
tial along a CCB direction helps to understand qualita- 
tively many of the properties of the tunnelling bounds. 
Suppose that the scalar fields take on values such that 



Ml 



- \ fi\ - \fj\ - 



(8) 



for one particular set of /, i, and j, where i,j = 1, 2, 3, 
are family indices with respect to the super-CKM basis. 
The other fields are all assumed to vanish. Then one 
can check that Vd = and the lowest potential energy 
subject to the above condition is 



V L . E . = [2 (Ml) 



f 



w 2 l « 2 
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(9) 



Due to the cubic term in a, a minimum can appear that 
is deeper than the SML vacuum unless 



U/,| 2 < A 



(10) 



By imposing this inequality, one can avoid a CCB mini- 
mum. 

This theoretical bound does not decouple even if 
sfermions grow heavier. On the contrary, it gets stronger. 
This property is most easily demonstrated in terms of a 
5 parameter. Using ([6]), one can translate (10) into the 
form 



with the superpotential 



U //./ /.,,V , e) + H d QiXfj d) - //„ (IX;, u) 
+ fiH u Hd- 



(5) 



The index a in pj) and ^ runs over each superfield in 
the model and <f> a is its scalar component. 



\(S{j)lr\ < m fn _ 



[2 (M 



2 )/ 
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(11) 



where rrif k is the mass of fermion This bound scales 
as the inverse power of the average sfermion mass. This 
behaviour is at variance with an FCNC limit on |(<^)l.r| 
that grows as sfermions become heavier. Another point 
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to notice is that the above restriction is independent of 
tan/3. As will be shown, this is not the case with a 
metastability bound. 

In addition to the CCB minima, UFB field directions 
can appear due to large trilinear couplings [3]. If one 
is worried only about the existence of a deeper point in 
the field space, both types of constraints should be taken 
into account. As far as a tunnelling process is concerned, 
however, a UFB direction is most likely irrelevant for the 
following reasons. As one traces this direction starting 
from the SML vacuum, the potential barrier is gcnerically 
thicker than along a CCB direction. Moreover, there is 
a section of path in which Vb does not vanish thereby 
making the barrier much higher as well. Therefore, a 
tunnelling process prefers a path close to a CCB direction 
that obeys Z?-flatness all the way. 

Even though expressed in a nice analytic form, the role 
of a CCB condition as a phenomenological constraint is 
rather unclear. Instead, the lifetime of the SML vacuum 
shall be required to be long enough. In a scmiclassical 
approximation |17j , the decay probability of a metastable 
vacuum per unit time per unit volume is given by 



T/V = A exp(-5[0]), 



(12) 



where A is a prefactor that needs to be guessed, S is the 
Euclidean action, and 4> is its bounce with 0(4) symme- 
try [18] . This symmetry allows one to use a single coor- 
dinate, i.e. the radial distance p, with respect to which 
the action can be written 



S[4>(p)] = 27T 2 



dpp 3 



dp 



V{cj>) 



(13) 



Note that the normalisation of the kinetic term is that 
of complex fields, as <j> is intended to be eventually the 
MSSM scalars. The field configuration 4>{p) is a station- 
ary point of S[(f>(p)] subject to the boundary conditions, 

^(p->oo) = ^+, 5r(p = o) = o, (14) 

dp 

where 4> + is the false vacuum. 

In the context of the MSSM, (f> represents the vector 
of scalar fields in the model, and 0+ is the metastable 
SML vacuum. One can evaluate the action by replacing 
V in ( p~3] ) with the scalar potential |l]) plus a constant 
term such that V (<(>+) = 0. 



On the other hand, the prefactor A in ( 12 ) is prac- 
tically impossible to compute. Based on dimensional 
grounds, it is guesstimated to be the fourth power of 
the mass scale appearing in the problem. In this work, 
A = (100 GeV) 4 shall be used with which one can esti- 
mate the smallest acceptable action. The lifetime of the 
whole observable universe staying at the SML vacuum 
is roughly V/(Tt^) where to ss 10 Gyr is the age of the 
universe. One can keep this longer than to by demanding 
that [TU] 



The above value of A may be regarded a little small in 
view of the sfermion masses in the TeV range. Choosing a 
small A implies a conservative parameter space exclusion. 
Nonetheless, the constraint is much less sensitive to the 
variation of A than of S. If desired, one could translate 
the multiplicative uncertainty of A into the additive one 
of S by taking the logarithm. For instance, taking A 
instead to be the fourth power of the reduced Planck 
scale would push up the lower bound on S by 150. 

Another type of uncertainty arises from radiative cor- 
rections. It is important to estimate their effects espe- 
cially because the numerical analysis uses the tree-level 
potential. One can implement the one-loop correction to 
the tunnelling rate simply by using the one-loop effective 
potential 



1 



VxL = V+-r-sVM+Sl, 



16tt 2 



(16) 



instead of the tree-level potential V in ( 13 ). The constant 
f2 should be adjusted so that Vih((j>+) = 0. The one- loop 
correction in the DR scheme is 1 19 



V™ =^(-l) 2s »(2 S „ + l) 



ln- 



ar) 



where m 2 is the field-dependent mass-squared eigenvalue 
of the n-th degree of freedom with spin s n . The renor- 
malisation scale Q should be chosen so that it minimises 
This means that the optimal scale is a function of 
<j>. Taking this scale would be satisfying but technically 
involved. For numerical computation, the scale shall be 
set equal to the common diagonal component of the soft 
sfermion masses. This is to be supplemented with error 
estimation. Fortunately, the bounce consists of <j> limited 
within a range not too far from the sfermion mass scale 
[7]. Therefore, V^ 1 ' does not diverge even if the CCB 
minimum is taken to infinity. 

One could estimate the uncertainty due to ignoring 
loop corrections to the action by evaluating the shift 
made by one-loop, 



Sih[<t> 



1LJ 



(18) 



where Sn, is given by ( |T3| with V replaced by the real 
part of the one- loop effective potential Vn in ( 16 ) [5D] , 



and the functional argument of each action is its bounce. 
One can approximate the difference ( 18 1 by 



AS = Sn[<P]-S[<j>] 



where 



2k* I dpp 3 Re[AU(</>(p))], 
/o 



AV = V 1L - V, 



(19) 



(20) 



S > 400. 



(15) 



keeping only the leading order term in AV. Note that 
AV is not simply equal to V^ /(16tt 2 ) + fl as may be 
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suggested by ( 16 ) since the parameters entering V and Therefore, negating /z can be compensated for by flipping 



Vil in (|20j) have different values. In particular, those pa- the sign of either Q 3 or u%. Note that this can be done 

without touching the trilinear term in question. In the 
case of a flavour-conserving trilinear, by contrast, the 
signs of A33 and \x cannot be disentangled [5T]. 

For faster computation, every scalar field is constrained 
to be real-valued. This is not completely general in that 
for a given set of real- valued {<j) a }, 



rameters pertinent to the electroweak symmetry breaking 
must deviate. Otherwise it would be impossible for both 
V and Vil to result in realistic Higgs vacuum expectation 
values. 



III. PROCEDURE OF ANALYSIS 

The method of numerical analysis is spelt out in detail. 
The same procedure is repeated for each of the species 
/ = e,d,u and for each of the six flavour-off-diagonal 
trilinear couplings. 

The families of the scalars / and f c are defined in 
the super-CKM basis where the corresponding Yukawa 
matrix Ajy is diagonal. The CKM mixing is neglected so 

that both A" and X d are diagonal at the same time. This 
simplifies the computation as the flavour-mixing Yukawa 
couplings are dropped from the superpotential. They are 
of 0(Vt s \t) ~ A 2 or less, and presumably are not very 
important in view of the fact that a coupling of similar 
order Af> ~ 0.1 does not make a sizeable difference as will 
be demonstrated later in this article. 

As the focus is on flavour- violation by the trilinear cou- 
plings, the soft masses shall be assumed to be in a uni- 
versal form, 



(21) 



Among the trilinear couplings, one A?- with i ^ j is 
scanned with the step size of 0.1 TeV while all the other 
components are kept zero. The diagonal, of no interest, 
is set to zero. One might well imagine a situation where 

more trilinears are turned on in addition to A/,. The 

l j 

additional couplings do not block an existing tunnelling 
path but may open more |21) . In this sense, the bounds 
obtained in this work should be regarded as generous. 

The phase of a single Afj can be rotated away by 
a lepton-flavour number transformation. Similarly, the 
phase of Afj or A^ may be absorbed into the squark fields 
by 'quark-family number' transformation in the limit of 
no CKM mixing that is taken here. 

The gauge and the Yukawa couplings are set at the 
renormalisation scale Q = my. Therefore, A^- presented 
later are also supposed to be at this scale. The other 
parameters shall be set as follows unless stated otherwise: 



tan/3 = 10, n = 0.5 TeV, m H a = 0.5 TeV, 

3 TeV. 



(22) 



The Higgs mass parameters, tit/jj , m? H , and 6, are deter- 
mined from these by the tree-level electroweak symme- 
try breaking condition. More parameter space is also ex- 
plored by varying each of the above. A remark is in order 
regarding the sign of ji. The only non- negligible term af- 
fected by this in the scalar potential is 2 Re[/iHdQlXtu"]. 



mm V(e 1 ' 
o<e a <2TT 



t>a) < 



min V(e l 

' =0,7T 



(23) 



even if all the parameters in V are real. One may worry 
that there can be a complex- valued path in the field space 
that costs less potential energy. Nevertheless, it can be 
checked that the equality holds as long as {4> a } stays 
sufficiently close to one of the directions in ^ , and thus 
the potential is dominated by a few terms including the 
activated trilinear term in particular. It turns out that 
the numerical bounce does reveal this property. 

Under the same condition, one can show that it costs 
less potential energy for the squark colour directions to 
be aligned. This justifies retaining only one colour com- 
ponent of each squark assuming that the others are all 
zero. 

As the problem is concerned with a flavour-off-diagonal 
coupling, one should put at least two families of scalar 
fields into action as variable degrees of freedom in ad- 
dition to the Higgs doublets. For the lepton sector, the 
following 10 variables are used: H^, H u , Li^e?, where the 
index i runs over the two families involved. Right-handed 
sneutrinos are not taken into account. The quark sector 
needs two more variables: i/^, H u ,Qi, d?, u\. 

The bounce is computed numerically. The Euclidean 
action (13 1 is approximated by a function of the fields 



that are put on the discretised p axis. Each lattice point 
has an independent set of the aforementioned 10 or 12 
degrees of freedom. One can find a stationary point con- 
strained by the boundary conditions ( 14 ) using a method 



derived from that in Ref [22 . The original series of steps 
had to be modified since it is not suitable for a prob- 
lem in which the true vacuum is very far from the false 
vacuum in comparison to the thickness of the potential 
barrier. This happens if the Yukawa couplings in ^ are 
small which is mostly the case except for X t . The details 
of the revision are presented in a separate article |23j . 
With this method, one does not have to know where the 
true minimum is, or even whether or not the potential is 
bounded from below. The price to pay is that one needs 
to minimise a function obeying constraints given in im- 
plicit forms. This can be done for instance by using the 
Ipopt package [24] . 

At the first stage of the computation, one needs to 
choose an initial field configuration. The above method 
needs a trajectory in the field space that connects <fi+ to 
any point 4> e over the barrier such that V{4> e ) = V{4>+). 
To find a <f> e , a point leaves the SML vacuum and hikes 
along a valley in the direction of increasing \<p\ until it 
goes down and touches the initial altitude. 
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TABLE I. The Euclidean action S decomposed into parts 
coming from the kinetic term, Vd, Vf, V so ft, and V, respec- 
tively labelled St, Sd, Sf, Ssoft, and Sv, in each of the three 
cases shown in Figs. [2] The last row shows the estimated 
one-loop effect. 



FIG. 1. Upper bounds on flavour-violating as functions 
of mj for / = e,d,u and i,j — 1,2,3. They correspond to 
the following choice of parameters: tan/3 = 10, fi = 0.5 TeV, 
m H o =0.5 TeV. Each group of four vertical lines shares the 
same value of mj given by the average of their horizontal 
coordinates. 



Finally, the SML vacuum lifetime is compared with 



the age of the universe using ( 15 1, to accept or reject the 
given parameter set. 

Apart from the scanning procedure described above, 
error estimation is performed for a few selected cases by 
evaluating the approximate one- loop effect ( 19 ). The pa- 
rameters entering the effective potential ( 16 ) are gener- 
ated by softsusy [35] with options set for one-loop Higgs 
masses and tadpolcsj^] Unlike the tree-level potential, 
Uil depends on the gaugino masses as well. They are 
picked up so that Mi. 2 , 3 = mj for order-of-magnitude 
estimation. 



IV. UPPER BOUNDS AND 
INTERPRETATIONS 



Requiring that the bounce action obey the condi- 
tion (15), one gets the upper bounds on the off-diagonal 
trilinear couplings shown in Fig. [T] It turns out that 
the six components of A\- have almost the same bounds 
(falling within 100 GeV from one another), and therefore 
a single vertical line is used to display them all together. 
The same is true of the six components of Afj. As for 
A™j, the upper bounds highly depend on whether either 
of the family indices is 3 or not. The bounds on A\ 2 and 
A 21 can be displayed together, and likewise the other four 
couplings. 

One can notice several features on the plot. First of 
all, the upper bound on an off-diagonal A grows as the 



the bound on Af, is 
This fam- 



corresponding sfermion mass is increased. The depen- 
dence turns out to be almost linear although not exactly. 
Second, the bound is symmetric under the interchange 
of the two family indices. Third 
independent of i and j as mentioned above 
ily independence is discovered also in Afj . This property 
contrasts with the Yukawa dependence of the stability 
bounds (10). Fourth, the bound on A\^ is nearly the 

Fifth, the bound on A\ 2 is slightly 

d 



same as that on Af 4 



1 There is a bug in softsusy version 3.1.7 that discards even one- 
loop tadpoles when the option mimRewsbLoops is set to 1. This 
problem does not appear with the default setting. 



lower than that on A^ . Sixth, the bound on A" 3 23 is 
higher than that on A" 2 and also tend to be higher than 
the other bounds for mj > 1 TeV. Interpretations of 
these observations shall be given in what follows. 

For understanding the results, it should be instruc- 
tive to examine contributions from different components 
of the scalar potential separately. Three instances are 
shown in Figs.[2j one with lepton flavour violation (LFV) 
and the other two with up-type quark flavour change. 
The figure for the down-quark sector is nearly identical 
to Fig. [2]ja) and therefore has been omitted. Each plot 
displays the differential action that arises from each of 
Vbj Vp, V so it, and their sum denoted by V, as one moves 
along the trajectory of the bounce profile. The horizontal 
axis is the distance traversed in the field space starting 
from the point of p = 0. Every curve arrives on the right 
at the false vacuum. Therefore, the area under (or over) 
each curve is the action coming from the given potential 
component. Also shown is the (positive-definite) differ- 
ential action from the kinetic term. Obviously, the total 
Euclidean action is given by the sum of the area under 
the light grey curve and the area of the light grey region 
(yellow online) . The parameters used in Figs. [2] are as 
in p2| besides the trilinear couplings indicated in the 
caption. Each parameter set leads to the action S ~ 400, 
that is, the A coupling is around the metastability limit. 

On each plot, the area under the Vr> curve is tiny com- 
pared to that under V. This means that Vd does not play 
a very important role in determining the action among 
the different components of the total scalar potential V. 
This point appears more concrete in Table [I] which lists 
the numerical value of each area. Examining the field 
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FIG. 2. Contributions of different potential components to the Euclidean bounce action. The nonzero trilinear coupling is (a) 
A23 = 4 TeV, (b) A12 = 3.8 TeV, or (c) Ah 1 -} = 5.3 TeV. The (signed) light grey areas (yellow online) display the action coming 
from the total potential energy density V that is the sum of Vr>, Vf, and V so f t , represented by the dashed (green online), the 
dot-dashed (blue online), and the solid grey (red online) curves, respectively. A black dotted curve is the estimated one-loop 
correction. The kinetic term is also displayed as a light grey curve. 



configuration, one can indeed find that the trajectory of 
the solution almost follows one of the D-flat directions 
in (J8|. The same holds in every case considered in this 
work in either the lepton or the quark sector. This al- 
lows one to understand the qualitative features of the 
metastability bound on a trilinear coupling by looking at 
the potential One property of Vl.e. is that it is sym- 
metric under the interchange of i and j if both (M% v ){j 
and A?j are symmetric. This makes it understandable 

that the bound on is symmetric. However, this ap- 
parent symmetry stems from the assumed form of the 
sfermion mass matrix in (21). Nevertheless, one could 



make a straightforward guess at the behaviour for a more 
general mass matrix: the bound should be a function of 
(M^ v )jy This property has been verified numerically. 

The second smallest contribution comes from Vf ■ It is 
totally negligible in a case for Af- or Af^. This explains 
why the upper bounds on Afj and Afj nearly coincide 
in Fig. [T] They do not depend on i or j, and are even 
independent of whether they are for charged sleptons or 
down-type squarks. The only differences between the 
two species are their gauge couplings provided that their 
masses and trilinear couplings are equal and that their 
Yukawa couplings are ignored. Since both Sf and So are 
negligible, the Yukawa and the gauge couplings do not 
make a noticeable difference between the bounds on A. 
and Afj. This should be clear from the potential in (|9 



with the Yukawa terms discarded. As for A" 2 , Sf is no 
as small but still minor compared to S so ft- Its bound in 
Fig. [l] is quite close to but a little stronger than those on 
A^ d , although the charm Yukawa is small as well. This 
split stems from the fact that A\ 2 is associated with the 
up-type Higgs unlike A\\ . This point shall be elaborated 



on later. 

A trilinear coupling including a stop is allowed to be 
substantially larger than A" 2 even though they involve 
the same Higgs field, as can be found in Fig. [T] The 
reason is the large top Yukawa coupling. It is evident 
in ([9| that Yukawas lift the scalar potential as quartic 
couplings, and so they slow down the tunnelling process 
if they are large. Therefore, large Xt implies weaker re- 
strictions on j4" 3 23. The nontrivial role of A* manifests 
itself in both Fig. pFc) and the rightmost column of Ta- 
bic |T] which show that Vp is not negligible as for the other 
A terms. 

The last element to be explained in Table |T] is the esti- 
mate of one-loop effects, which was defined in ( |19[ ). The 
differential version is plotted in Figs. [2| The magnitude 
of AS turns out to be roughly 3% of the tree-level S in 
the case of A 23 and is less significant in the other two. By 
full computation of the one- loop bounce 1L , the validity 
of AS as an approximation of the one-loop shift in ( 18 ) 



has also been checked. The latter quantity is indeed of 
the same order as AS for any case shown in the table. 
In combination with the smallness of AS, this should a 
posteriori support the reliability of an analysis employing 
the tree-level effective potential. 

Nevertheless, one should like to translate the variation 
in S into that in the maximum Afj. In any case, there are 
non- vanishing uncertainties coming from loop corrections 
as well as from the ignorance of the prefactor. Also, 
one could choose to impose a lower limit on S slightly 
different from that in (15). All these variances could be 
encoded as a constant added to S. Then, its influence 
on a trilinear coupling can be found in Figs. [3] If one 
allows for a change in S by a factor between 0.6 and 1.6 
for instance, the bound on A[- is altered by ±20%' anc ^ 
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FIG. 3. The plots show how the upper bounds change depending on the condition imposed on S. On each curve, the value of 
A{j leads to the corresponding S shown on the right, where / = e,d,u and i,j = 1,2,3. The following values of parameters 
were used: tan/3 = 10, /j, = 0.5 TeV, m H o = 0.5 TeV. 



each metastability bound presented later is weakened or 
strengthened by the same fraction. 

As is clear in Fig. [lj the tunnelling process gets slower 
as the squark and the slepton masses increase since they 
lift the potential ([9| quadratically. By the same token, 
one may envisage a similar effect from higher masses of 
Higgs fields, as either of them also acquires a vacuum ex- 
pectation value along the CCB direction Concern- 
ing the Higgs bosons, the parameters directly linked to 
the tunnelling are not the soft masses appearing in ^ 
(which can be tachyonic) but the physical mass eigen- 
values which describe the positive curvature of the po- 
tential around the SML vacuum. In the MSSM, one can 
specify the tree-level Higgs potential by choosing tan/3, 
77i //o , and \i. Once the former two parameters fix the 
complete Higgs mass spectrum, varying fj, makes no dif- 
ference therein. The dependence on each of these three 
parameters shall be discussed below. 

First, and the sfermion masses are varied while 

the other parameters are fixed according to ( 22 ) . The re- 



sults are set out in Figs.[4]as grey contours. As expected, 
enlarging m^a does provide more room for the trilincar 
couplings. However, the effect is significant only on those 
involving the down- type Higgs, displayed in Fig. [4][a). 
Contours of maximum Afj in Figs. Qb) and (c), do not 
exhibit a noticeable dependence on ttt^o. One can un- 
derstand this by considering Higgs mixing. The real part 
of each Higgs field can be expressed as (see e.g. Ref. [21)] ) 

Re#° « v u + (sin/3 h° - cos/3 H°)/V2, (24a) 
ReH% v d + (sin/3 H° + cos/3 h°)/V2, (24b) 

in terms of the mass eigenstates. The Higgs mixing angle 
a has been approximated by /3 — tt/2 since m H o 3> rriz 
in most of the parameter space considered in this work. 



For tan/3 = 10, H° and H% are mostly h° and H°, re- 
spectively. Since the mass of the lighter CP-even neutral 
Higgs boson is always restricted around mz, the potential 
curvature along the H® direction stays essentially con- 
stant in comparison to the large variation in the squark 
mass. This property should remain true at least qualita- 
tively for any tan/3 bigger than a few. 

A metastability bound should be easier to use if one 
has its algebraic expression. Within the parameter region 
shown in the figures, the bounds can be approximated by 
the empirical inequalities, 



Alf < 1.61m j 



0.59 m H o 

0.93„,0.07 



0.45 mf^m^l' + 37 GeV, 



Af 221 < 1.27 rrig- 44 GeV, 



32 GeV. 



^13,31,23,32 < 1-76 m g 



(25a) 

(25b) 
(25c) 



These are the least-squares fits to the numerical data 
,■ " i ! ' i i i 1 1 ( - 1 ! ' 1 1 : j i d. They are 



using the form, amr+ bra jjo +cm p ~ m l p 



f 



plotted as dashed black curves in Figs.|4]which reveal that 
the fits are indeed faithful approximations of the data. 
These inequalities should remain nearly the same even 
for different values of tan/3 and /i since the bounds are 
insensitive to them as will be shown below. It should also 
be remembered that m~ should be replaced by (M^ v ){ 

if the diagonal soft masses are highly nondegenerate. 

Next, the dependence on tan/3 is illustrated in Figs. [5] 
There are two origins of change as tan/3 varies. One is 
the Higgs mixing mentioned above. This works in the 
same way on all the species tha t cou ple to each type of 
Higgs. As tan/3 increases, H® in (24b) acquires more frac- 



tion of the heavier Higgs which impedes the tunnelling 
caused by A^ d . This effect appears as the slight rise of 
the points for higher tan/3 in Figs.^a) and (b). For 
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FIG. 4. Contours of upper bounds on flavour-violating on the (rnj,m H o) plane for / = e,d,u and i,j = 1,2,3. The 
grey contours [which look wiggly in Fig. (a)] are based on the data from grid scanning and the dashed black curves are the 
least-squares fits. They correspond to the following choice of parameters: tan/3 = 10, /i = 0.5 TeV. 
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3.5 



A/1 13,31,23,32 
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b 4.5 



4.0 
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10 
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FIG. 5. Upper bounds on flavour-violating A{- as functions of tan/3 for / = e,d,u and i, j = 1,2,3. They correspond to the 
following choice of parameters: /i — 0.5 TeV, m H o = 0.5 TeV, mj = 3 TeV. Each pair of vertical lines shares the same tan/3. 



Afp the change is in the opposite way but too small to 30, which does not seem very significant, 
be visible on the points for A\ 2 in Fig. [5jc). Remember 
that increasing m#o renders the tunnelling less probable. 
Therefore, the influence of Higgs mixing should be en- 
hanced by higher m^o . The other origin is the variation 
of Yukawa couplings, which dominates the tan/3 depen- 
dence of the trilinear couplings involving the third family. 
A higher value of tan/3 implies larger down-type quark 
and charged lepton Yukawas. This gives rise to slower 
tunnelling as explained before. This splits to some ex- 
tent the upper bounds that were close together for lower 
tan/3. The up-type Yukawas do not change as much as 
the down- type ones, and this leads to the milder depen- 
dence of A^ 3 23 in Fig.lBTc). After all, the alteration is at 
most 10%, found in A^, as tan/3 is increased from 3 to 



Finally, one should discuss the dependency on /i. As 
already stated, the Higgs mass eigenvalues and mixing 
angle do not rely on [i after tan/3 and m#o have been 
fixed. Therefore, the metastability bounds should remain 
the same while /j, is altered. This has been checked by 
trying two values of fi of different orders of magnitude: 
0.5 TeV and 3.0 TeV. For each /x, rnjyo has been scanned 
from 0.5 TeV to 5.0 TeV, with the other parameters set 



as (22). At these parameter points, the two choices of /j, 
do 



ead to identical upper bounds. 
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Mode 

fl — > 67 
t — y e*y 
t — > flj 



Current bound Future bound 



mates in Ref. [5]. An estimate in this reference is in the 



1.2 • 1CT 11 [27] 

3.3 • 10" 8 [29] 

4.4 • 10" 8 \M 



10" 
10" 
10" 



28 

m 
m 



TABLE II. Present and future experimental sensitivities to 
the branching ratios of the lepton flavour violating decay 
modes that probe (S^lr- 



MI 


Current bound 


Future bound 


(3i2)lR,RL 


1.4- 10" 3 [32] 




{8i 3 )hR,RL 


2.9- 10" 2 [32133] 


2- 10" 3 [5] 


{82z)lr,rl 


1.4- 10" 2 [32] 


5- 10 -3 [5] 


{8i 2 )lr,rl 


1.6 ■ 10~ 2 E2 





TABLE III. Experimental sensitivities to flavour-violating 
mass insertions in the squark sector for = Mz = 1 TeV. 
Expectations of future searches are quoted if available. 



V. IMPLICATIONS FOR FLAVOUR PHYSICS 

Having obtained the limits on flavour-violating trilin- 
ear couplings, one should consider what consequences 
they have for physical processes. In what follows, the 
metastability bounds shall be contrasted with other exist- 
ing constraints. They are based primarily on new physics 
searches within FCNC processes, of which prospects at 
running and planned experiments are also discussed. 
Other indirect bounds are compared together such as 
those coming from the vacuum stability, the p parameter, 
and naturalness of the CKM matrix. In this section, the 
degree of flavour- violation is expressed as a S parameter, 
defined in (J6J), that is normally used for flavour physics. 

The strongest experimental constraints on the LR 
mass insertions in the slepton sector come from radiative 
LFV decays. The restrictions on their branching ratios 
are collected in Table |H] The current upper limits are all 
at 90% confidence level (CL). The quoted future sensi- 
tivity to n — > ej is anticipated at 90% CL from the MEG 
experiment [28] . According to Ref. [5j , a SFF is expected 
to provide the displayed 90% CL bounds on t — > cy and 
r — > /17, whereas Ref. [3U] projects the r — > pj bound of 
3-10~ 9 . 

Each bound is translated into that on the correspond- 
ing mass insertion and presented below. For this, one 
should choose the other parameters on which the LFV 
decay rate depends. In the case of (Sfj)LR or (5^)rl with 
j > i, fixing the bino mass Mi is enough to determine 
B(lj — > ?j7) in the mass insertion approximation [31]. As 
a demonstration, it shall be assumed that M\ = mj. 

In the squark sector, there are multiple observables 
that restrict a single mass insertion. These constraints 
are combined to give the 95% CL upper bound on the 
modulus of each mass insertion in Table |III[ Also pre- 
sented are the moduli that in the future can be recon- 
structed at the level of 3<r. They are based on the esti- 



form of a range and the values in Table III are those at 
the lower ends, i.e. the most optimistic ones. The limit 
depends on the squark and the gluino masses as it arises 
from gluino-squark loops. The numerical values shown 
in the table are for rriq = M3 = 1 TeV, and they are pro- 
portional to the squark mass as long as rriq/M 3 is fixed. 
This ratio shall be maintained in the following plots. 

The FCNC constraints presented below should be re- 
garded as rough estimates. They rely on parameters that 
do not affect a metastability bound. Already at one- 
loop level, they depend on gaugino masses which have 
been chosen somewhat arbitrarily as stated above. More- 
over, chirally-enhanced higher order corrections can sig- 
nificantly strengthen or weaken a flavour constraint [34] . 
These corrections shall be ignored that depend on yet 
more parameters. 

Given the aforementioned experimental information, 
one could start the comparative analysis by examining 
the charged leptons. In Figs. [6j the vacuum bounds are 
displayed along with the present and the projected sensi- 
tivities of FCNC searches. The vacuum stability border 
is given by (111, above which the corresponding trilin- 
ear term gives rise to tunnelling. If one abandons the 
requirement that the SML vacuum be absolutely stable, 
then the limit is relaxed up to the line above which the 
false vacuum decays too fast in comparison to the age of 
the universe. As the sfermion masses increase, both the 
stability and the metastability bounds become stronger. 
However, the latter depends on tan/3 unlike the former. 
As Fig. pfa) shows, maximum Afj allowed by the vac- 
uum lifetime is insensitive to tan/3. Therefore, ^ im- 
plies a tighter restriction on (5fj)LR for higher tan/3. In 
contrast, the FCNC constraints loosen as the sfermions 
grow heavier. 

More specifically, Fig. [6^a) is concerned with the fi — > e 
transition. Within the displayed range, the tunnelling 
process does not play a significant role in comparison to 
the fi — > 67 decay. The MEGA experiment has ruled out 
more region than metastability so long as mr < 20 TeV 
for tan/3 = 30. There is even no probability of tunnelling 
for up to mr fa 3 TeV, beyond which the SML vacuum 
lifetime is still permitted to be either finite or infinite. 
Allowing for an arbitrary mass insertion consistent with 
each vacuum bound, one can estimate what slepton mass 
range may be accessible through fi — > ej. The stability 
condition sets the higher end around 10 TeV. If one opens 
up the possibility of a long-lived false vacuum, the terri- 
tory extends up to mr~ 60 TeV where the metastability 
bound for tan/3 = 30 crosses the MEG sensitivity. As 
is clear in the plot, this crossing point depends on tan/3, 
and is about 100, 200 TeV for tan/3 = 10, 3, respectively. 
Above this mass scale, the mass insertion is suppressed 
too much to be observed. The metastability bound has 
been numerically computed only up to m~ = 5 TeV, and 
the preceding estimates are based on linear extrapolation. 

In a similar way, one can interpret Fig. |6jb) which 
shows the 1-3 and the 2-3 mixings at the same time. 
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FIG. 6. Comparison of constraints on (<5y)iij. The thick black lines with negative slopes are bounds from the requirement 
that S > 400, for three values of tan/3 represented by different dashing styles. The thick grey lines are the stability bounds. 
The thick dark grey lines (blue online) with positive slopes are the current bounds from flavour- changing processes. The thin 
dotted lines (blue online) are anticipated experimental bounds. 



Obviously, they are associated with r — > and r — > 
fj/y, respectively. In contrast to ji — > ey, their current 
experimental bounds are less powerful and leave ample 
room for vacuum metastability even for slepton masses 
as low as 500 GeV. If one were to insist on permanent 
stability, it could be an explanation for why LFV in r 
decays has not been discovered yet. This requirement 
would rule out even the possibility of observing r — > e~f 
or t — > /X7 at a SFF, provided that they arise solely from 
the trilinear terms and that mj > 600 GeV. Adopting 
the vacuum longevity constraint instead makes the LFV 
search more promising. A SFF might be able to find 
the virtual effect of a LR insertion involving the third 
family if sleptons are lighter than 1.0,1.6,2.7 TeV for 
tan/3 = 30, 10, 3, respectively. If one wishes to use the 
projected r — > ^7 bound of 3 • 10 -9 from Ref. [30], one 
may shift up the corresponding line in the plot by the 
factor of a/3/2. 

Next, the down- type quark sector is investigated. The 
structure of each of Figs. [7] is the same as a leptonic one. 
The left panel is devoted to the s — > d transition. The 
kaon physics constraint is weaker than the stability con- 
dition but stronger than metastability up to the mass 
scale around a few TeV. Again, the crossing position 
varies according to tan/3 which can be found in Fig. (7^a). 
If squarks are heavier than this mass scale, (Sf 2 )LR,RL 
are bounded to be smaller than the sensitivity of present 
measurements. However, one should keep in mind that 
kaon dynamics is plagued by large theoretical uncertain- 
ties and so is the experimental limit. 

The mass insertions relevant to the bottom quark are 
shown in Fig. [7](b) . The features of the 1-3 and the 2- 
3 sectors are similar, with differences arising from the 
FCNC constraints. The current B physics data does not 
give much information with respect to whether or not 



a vacuum decay can be triggered by the related trilin- 
ear term. As long as one stays within the stable region, 
the scope of a SFF stops around 1.4 TeV via the 1-3 
mixing and a little lower via 2-3. If one goes beyond 
the stability limit, the reach of new physics search in 
b — > d transitions is raised up to 2.1,3.4,5.9 TeV for 
tan/3 = 30, 10, 3, respectively. The respective b — > s 
counterparts are 1.3,2.2,3.7 TeV. As before, these lim- 
its are set by the age of the universe. Remember that 
there are hadronic uncertainties in heavy flavour physics 
as well albeit less than in the kaon sector. Also, the future 
bounds in Table [TTT| have rather big errors. In particular, 
new physics is much easier to unveil if there is an ex- 
tra complex phase, since CP violation is less affected by 
hadronic uncertainties in general. These points should 
be taken into account when one interprets the above nu- 
merical results. 

Before moving on to the up sector, it should be ap- 
propriate to recall the Higgs mass dependence. The pre- 
ceding graphs correspond to the choice mjia = 0.5 TeV. 
As demonstrated in Fig. |4]ja), higher ?ti#o makes bigger 
i&if) lr acceptable. This can extend the boundary of 
reachable mass upwards. 

Finally, Figs. [8] are allocated to flavour violation of the 
up-type squarks. One finds an outstanding distinction 
between the metastability bounds shown here and in the 
preceding plots, i.e. the way they depend on tan/3. This 
is because the up-type squarks couple to H u unlike the 
fields considered above. According to (J6J, the bound be- 
comes looser as tan/3 increases although the dependence 
is much milder. 

Regarding the 1-2 sector, there is a constraint from a 
flavour-changing process. In Fig.^a), the D°-D° mixing 
bound is plotted. It does not forbid the SML vacuum 
having a finite lifetime, but does restrict the tunnelling 
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FIG. 7. Comparison of constraints on (SfAhR- The thick black lines with negative slopes are bounds from the requirement 
that S > 400, for three values of tan/3 represented by different dashing styles. The thick grey lines are the stability bounds. 
The thick dark grey lines (blue online) with positive slopes are the current bounds from flavour- changing processes. The thin 
dotted lines (blue online) are anticipated experimental bounds. 




FIG. 8. Comparison of constraints on {SfAhR- The thick black lines are bounds from the requirement that S > 400, for three 
values of tan/3 represented by different dashing styles. The thick grey lines are the stability bounds. The thick dark grey line 
(blue online) in panel (a) is the current bound from D°-D° mixing. In panel (b), the thick dark grey lines (blue online) are 
the bounds from the p parameter. 



rate below the acceptable level for < 4 TeV. Beyond 
this mass scale, the metastability bound is stricter than 
the current data. 

There is not much experimental information on the 
LR insertions related to the top squark. In particular, 
no FCNC constraint is available yet. Instead, Fig. [8jb) 
employs another indirect limit coming from electroweak 
precision observables. It is obtained by requiring that the 
squark contribution to the p parameter [35] do not exceed 
5.5- 10 -4 [36]. It is based on the assumption that {SfytR 
is not correlated with (6% ) lr so that a large value of the 
latter breaks the custodial symmetry between the up- 
and down-type squarks significantly. As squarks grow 



heavier, this limit does not decouple either but tightens. 
However, the pace is slower than those of the vacuum 
constraints. As a result, the vacuum longevity bound is 
more stringent in most of the parameter space. Another 
point to notice is that the stability and the metastabil- 
ity curves are very close to each other compared to the 
foregoing cases that do not involve the top quark. This 
is mainly because the large top Yukawa coupling renders 
the inequality (10) less restrictive. In the plot is shown 



even a small region which is ruled out by the tunnelling 
constraint for some tan/3 even though it satisfies the sta- 
bility requirement. Of course, a classically stable vacuum 
cannot decay This apparent nonsense stems from the 
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MI 


Metastability 


Naturalness 




0.026 


0.0011 


{8*3) LR 


0.026 


0.0010 


(823) LR 


0.026 


0.010 


{812) LR 


0.21 


0.011 


(813) LR 


0.29 


0.062 


(823) LR 


0.29 


0.59 



TABLE IV. Metastability and naturalness bounds for mj = 
M3 — 1 TeV. The metastability bounds on (Sf^LR are pro- 
portional to cos/3 and the displayed values are for tan/3 = 10. 



fact that the condition (10), used for drawing the curve, 
is not optimal [21] , That is, obeying this condition is 
not enough to exclude any CCB minimum to which tun- 
nelling can occur. It is beyond the scope of this article 
to pin down the necessary and sufficient conditions for 
vacuum stability with respect to the flavour-violating tri- 
linear terms. In any case, the stability bound is not the 
primary concern here. 

There is another class of theoretical bounds that scale 
as the inverse power of the sparticle masses. They 
are naturalness bounds that act on half of the flavour- 
changing squark mass insertions, (£y U )z,ii with i < j 
|37j . They are based on the requirement that the super- 
symmetric loop corrections do not exceed the measured 
values of the CKM matrix elements. They are compared 
with the tunnelling limits in Table |IV| One can notice 
that metastability is tighter on (6% 3 )lr an d that the two 
types of constraints are comparable on (S^lr especially 
for high tan/3. For the remaining four insertions, nat- 
uralness tends to be stronger. This strength depends 
on the degree of fine-tuning that one is willing to allow, 
which by contrast has nothing to do with the metastabil- 
ity bounds. Note that the slepton sector and the other 
half of the squark sector LR insertions are not restricted 
by the naturalness in renormalisation of the lepton or 
the quark mixing matrix. These insertions can be con- 
strained by two-loop corrections to light fermion masses 
but only in a combination with another insertion |38j . 
These multiple insertion bounds shall be left out of the 
present comparison as the vacuum lifetime limits a LR 



insertion even if it is the only non- vanishing component. 



VI. CONCLUSIONS 

In the context of the MSSM, upper bounds on the 
flavour-violating trilinear soft terms have been obtained 
by demanding that the standard vacuum be long-lived. 
Obviously, metastability leaves more room than stability. 
As with the latter, however, the former results in a limit 
that does not decouple even if raj increases. A distinct 
property of the new bounds is that those on the 12 com- 
ponents A^j d are nearly the same. This is because the 
bounce almost follows a D-flat direction and the corre- 
sponding Yukawa couplings are negligible. Dependence 
on other parameters has been reported as well such as 
Higgs masses, tan/3, and li. The bounds are fairly sta- 
ble against radiative corrections in view of the estimated 
errors. 

Prospects for indirect new physics discovery have been 
discussed in a scenario where flavour-violating trilinear 
terms can be arbitrarily large within the limits from vac- 
uum longevity. Being sensitive to mass ranges as high 
as 60 TeV or more, \x — > could be the most powerful 
probe of new physics. The potential of a SFF has also 
been considered. It should be able to cover up to a few 
TeV through B physics and LFV r decays. Lower tan/3 
tends to extend the territory, but it does not seem to grow 
an order of magnitude beyond the reach of the LHC. It 
should be reminded that these conclusions are based on 
the assumption that the LR insertions are the only source 
of extra flavour violation. As for the stop couplings, 
metastability provides bounds that are stronger than any 
existing experimental constraints even for rather light 
squarks. 
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